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ELECTROMAGNETIC  SHIELDING  EFFECTIVENESS  FOR 
ISOTROPIC  AND  ANISOTROPIC  MATERIALS 


I.  INTRODUCTION 

The  purpose  of  this  report  is  to  present  analysis  techniques  applicable 
to  the  assessment  of  electromagnetic  shielding  effectiveness  of  a  variety  of 
types  of  materials  including  advanced  composite  materials  (e.g.,  fiber-rein¬ 
forced  epoxies  where  the  fibers  are  "long"  and  "oriented"),  conductive  filled 
thermoplastics  (randomly  oriented,  chopped  fibers  or  other  conductive  particles) 
and  the  usual  metallic  shields.  Advanced  composite  materials  are  now  widely 
used  as  structural  and  surface  components  in  aircraft  and  spacecraft.  Chopped- 
fiber-filled  structural  foams  are  used  as  enclosures  for  electronic  equipment. 
Such  applications  require  accurate  assessment  of  electromagnetic  shielding 
effectiveness . 

The  following  sections  present  a  unified  approach  to  plane-wave  shield¬ 
ing  analysis  using  ABCD  and  scattering  parameter  techniques  based  on  analogies 
with  distributed  networks.  Related  techniques  and  preliminary  results  were 
presented  in  an  earlier  report  [1].  Some  key  results  from  [1]  are  included  in 
this  report  for  completeness. 

Chapters  2  and  3  present  background  material.  Several  important  results 
for  isotropic  materials  are  developed  in  Chapter  3  and  presented  in  new  formats 
Chapter  4  introduces  scattering  parameter  techniques.  Chapter  5  treats  aniso¬ 
tropic  multilayered  shields  from  several  points  of  view,  including  a  very 
general  arbitrary  polarization,  arbitrary  angle  of  incidence  analysis. 


2 .  BACKGROUND 

An  electromagnetic  shield  is  an  enclosure  intended  to  prevent  electro¬ 
magnetic  energy  within  it  from  escaping  and  conversely  to  prevent  electromag¬ 
netic  energy  exterior  to  it  from  penetrating  to  the  interior.  A  perfect  shield 
would  provide  total  isolation  between  interior  and  exterior.  In  practice  the 
degree  of  isolation  between  interior  and  exterior  is  less  than  perfect  and 
seldom  exceeds  100  dB. 

The  degret  of  isolation  provided  between  interior  and  exterior  of  a 
shield  is  commonly  called  its  shielding  effectiveness.  A  measure  of  shielding 
effectiveness  could  be  obtained  by  determining  field  strength  at  a  particular 
location  inside  the  shield  and  then  removing  the  shield  and  determining  the 
ambient  field  strength.  Errors  can  occur  in  such  measurements  due  to  such 
effects  as  modification  of  the  sources  of  the  external  ambient  field  by  the 
presence  of  the  shield,  and  proximity  effects  on  the  internal  field  measuring 
device  caused  by  the  shield  walls. 

Shielding  effectiveness  is,  thus,  the  insertion  loss  incurred  by  elec¬ 
tromagnetic  energy  in  passing  from  an  input  medium  through  a  shield  into  an 
output  medium.  Both  reflection  and  absorption  contribute  to  shielding  effec¬ 
tiveness.  Aircraft  require  substantial  electromagnetic  shielding  to  protect 
sensitive  internal  electronics  from  extraneous  signals.  Some  vital  areas 
particularly  sensitive  to  inadequate  shielding  include  low-level  integrated 


circuit  communication  and  navigation  equipment,  fly-by-wire  systems,  fire 
control  systems,  and  electro-explosive  devices. 

In  a  typical  situation  the  skin  of  an  aircraft  might  be  considered  the 
shield  with  the  air  external  to  the  vehicle  as  the  input  medium  and  a  complex 
mixture  of  personnel,  air,  cables,  electronic  devices,  fittings,  etc.  of  the 
aircraft’s  interior  as  the  output  region.  The  external  electromagnetic  field 
to  be  excluded  from  the  interior  of  the  shield  might  be  produced  by  any  of  a 
variety  of  electromagnetic  sources  including  a  direct  lightning  strike  to  the 
aircraft,  a  nearby  lightning  strike,  a  nearby  radar  or  other  high-powered 
transmitter,  or  perhaps  by  a  nuclear  detonation  (EMP) .  Alternatively,  the 
shield  may  be  the  cabinet  enclosure  of  an  on-board  piece  of  equipment  whose 
interior  is  to  be  protected  from  all  the  above  types  of  signals  plus  possible 
interference  from  other  on-board  equipment. 

Potential  interfering  electromagnetic  fields  may  arise  from  "nearby"  or 
"distant"  generating  sources.  "Near"  to  a  generating  source,  usually  either 
the  electric  or  the  magnetic  component  of  the  field  will  dominate.  For  sources 
that  resemble  loops  of  current,  magnetic  field  dominates  in  the  region  near  to 
the  source.  For  sources  that  resemble  linear  dipoles,  electric  field  dominates 
in  the  region  near  to  the  source.  Sufficiently  far  from  either  type  of  source, 
the  propagating  field  becomes  a  plane  wave  in  which  the  field  energy  is  equally 
divided  between  the  magnetic  and  electric  field  components. 

Wave  impedance  is  defined  as  a  ratio  of  electric  to  magnetic  field 
components.  For  a  plane  wave  that  ratio  is  120tt  ohms  or  approximately  377  ohms. 
Near  to  current  loop  type  sources  the  wave  impedance  is  very  small  compared  to 
377  ohms  because  the  magnetic  field  is  large  relative  to  the  electric  field. 
Near  to  linear  dipole  type  sources  the  wave  impedance  is  very  high  compared  to 
377  ohms  because  the  electric  field  is  large  relative  to  the  magnetic  field. 
Wave  impedance  of  the  energy  incident  on  a  shield  is  a  critical  factor  in 
determining  whether  or  not  the  shield  will  be  effective  in  excluding  the  signal. 
It  turns  out  (as  will  be  seen)  that  it  is  much  more  difficult  to  shield  "low- 
impedance"  magnet ic-type  waves  than  to  shield  either  "high-impedance"  electric- 
type  waves  or  "normal-impedance"  plane  waves.  Shielding  against  low- impedance 
waves  is  particularly  difficult  at  low  frequencies. 

A  rough  guideline  for  separating  sources  into  "nearby"  and  "distant" 
types  is  to  consider  it  a  nearby  source  if  it  is  closer  than  one-tenth  of  a 
wavelength  to  the  shield.  At  100  KHz,  1  MHz,  100  MHz,  and  1  GHz  one- tenth 
wavelength  is  approximately  1000  feet,  100  feet,  1  foot  and  0.1  feet,  respec¬ 
tively.  Thus,  for  frequencies  up  to  about  1  GHz,  one-tenth  wavelength  is  a 
relatively  large  distance  compared  to  the  size  of  modern  circuits  and  "near 
fields,"  at  least  from  other  on-board  equipment,  may  cut  across  several 
circuits . 

Practical  shields  seldom  have  shielding  effectivenss  greater  than  100 
dB.  Metals  are  in  general  good  electromagnetic  shields.  Shielding  effective¬ 
ness  of  a  metal  structure  is  degraded  by  the  presence  of  fabrication  seams  and 
joints,  access  doors  and  windows  and  other  apertures.  A  metal  aircraft  with 
the  usual  apertures  seams  and  cracks  will  have  an  effective  overall  shielding 
of  the  order  of  20  dB  at  UHF  frequencies.  Commonly  used  composite  materials. 
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composed  typically  of  small,  relatively  poorly  conducting  fibers  embedded  in  an 
insulating  matrix,  provide  considerably  less  shielding  particularly  at  lower 
frequencies.  Metal  matrix  composites  and  insulating  matrix  composites  with 
metallic  fibers  are  being  studied,  and  from  the  point  of  view  of  electromag¬ 
netic  shielding  such  materials  would  clearly  be  superior. 


3.  SHIELDING  THEORY 

Shielding  effectiveness  (S.E.)  is  an  exceptionally  difficult  quantity  to 
evaluate  for  a  given  material  because  S.E,  dpends  not  only  on  intrinsic  material 
parameters  but  also  is  a  strong  function  of  shield  geometry  and  of  both  inter¬ 
nal  and  external  environments  of  the  shield.  In  actual  situations  internal 
and  external  shield  environments  are  usually  complicated  and  time  varying  (e.g,, 
internal-personnel  and  equipment  inside  an  aircraft;  external-aircraft  in 
hanger,  on  runway,  in  air).  In  addition  real  shields  have  a  variety  of  seams, 
joints  and  apertures  that  generally  degrade  shielding  effectiveness.  Approxi¬ 
mations  are  necessary  in  order  to  reduce  a  shielding  configuration  to  a 
manageable  electromagnetic  boundary  value  problem. 

Plane  wave  shielding  theory  has  long  been  used  as  an  aid  in  character¬ 
izing  materials  and  providing  a  baseline  measure  of  shielding  effectiveness. 

In  this  report  plane  wave  shielding  theory  is  extended  to  several  classes  of 
anisotropic  materials  and  recast  in  the  very  useful  scattering  parameter  format . 
Coupling  mechanisms  that  degrade  instrinsic  shielding  are  briefly  described. 

3 . 1  Plane  Wave  Normally  Incident  on  Infinite  Flat  Plate 

The  manner  in  which  an  electromagnetic  shield  transmits  plane  electro¬ 
magnetic  waves  has  been  shown  [2]  to  be  analogous  to  the  manner  in  which  a 
conventional  transmission  line  transmits  electrical  current  and  voltage.  An 
idealized  configuration  consisting  of  a  plane  wave  normally  incident  on  an 
infinite  flat  plate  shield,  as  shown  in  Figure  1,  provides  a  useful  baseline  or 
reference  value  of  shielding  effectiveness  <>  It  also  serves  as  a  surprisingly 
good  approximate  model  for  a  number  of  more  complex  configurations  that  are  used 
for  laboratory  evaluation  of  shielding  effectiveness.  As  a  first  step  it  is 
assumed  that  the  shield  material  is  homogeneous  and  is  isotropic  in  the  plane 
of  the  shield.  This  hypothesis  obviously  includes  the  usual  metallic  shields 
but  it  also  appears  to  include  multilayer  graphite  epoxy  laminates  in  which 
the  layers  are  oriented  at  different  angles  (e.g.,  0°-90°-0°).  Unidirectional 
laminates  and  perhaps  even  multilayer,  mixed  orientation  boron  laminates  require 
a  different  hypothesis  taking  into  account  the  three-dimensional  anisotropy 
effects.  The  less  complex  behavior  of  the  multilayer,  mixed-orientation 
graphite  laminates  arises  from  the  f iber-to-f iber  contact  between  layers.  The 
more  complex  unidirectional  layer  configuration  has  also  been  modelled  and  is 
discussed  in  Section  5, 3  of  this  report. 

Using  the  well  known  transmission-line  analogy  [2]  of  plane  wave  propa¬ 
gation  as  shown  schematically  in  Figure  2,  all  types  of  2-port  representations 
[31  may  be  utilized  to  represent  the  shield.  The  "Zn  and  MABCDM  parameter 
representations  have  been  widely  and  successfully  used  in  the  past.  Scattering 
("S")  parameters  [4]  offer  in  many  instances  distinct  advantages  such  as  ease 
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of  automated  measurement  and  direct  interpretation  in  terms  of  transmitted  and 
reflected  components.  Scattering  parameters  are  treated  in  Chapter  4  and 
applied  to  multilayer  anisotropic  shields  in  Chapter  5.  The  scattering  param¬ 
eter  approach  is  equally  valid  and  easy  to  use  for  the  isotropic,  homogeneous 
shield.  However,  to  facilitate  comparison  with  certain  classical  results,  the 
immediately  following  analysis  is  carried  out  in  terms  of  ABCD  parameters.  A 
simple  scheme  for  conversion  to  S-parameters  is  given  in  Section  4.3.  The 
computer  programs  used  for  numerical  evaluation  of  shielding  effectiveness  use 
a  combination  of  multiport  ABCD  and  scattering  parameters. 

Figure  3  (a-c)  illustrates  the  details  of  the  schematic  representation 
of  plane  wave  propagation  via  the  transmission-line  analogy.  Proper  choice  of 
the  source  and  load  wave  impedances  (Z  and  Z  )  permit  this  model  to  be  used 
even  when  the  incident  wave  is  not  a  pYane  wave,  e.g.,  in  the  near  field  of  a 
loop  of  current.  This  technique  is  illustrated  for  a  two-loop  measurement 
configuration  in  Section  3.2.  The  basic  procedure  is  to  analyze  the  trans¬ 
mission  line  model  to  determine  its  total  insertion  loss.  That  insertion  loss 
is  then  the  same  as  the  shielding  effectiveness  of  the  shield  modelled. 

The  overall  ABCD  matrix  for  the  shield  with  the  effects  of  the  output 
medium  included  is  obtained  as  the  product  of  the  matrices  for  the  individual 
components  and  is  given  by: 


where  A  =  D  =  coshO  and 

a, 

B  =  H  sinh0 
C  =  1/n  sinh0 


n  =  Ai/e  =  intrinsic  wave  impedance 
0  =  yt  =  electrical  length 
Y  =  propagation  factor 
y  =  permeability 
£  =  permittivity. 


Thus , 
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Figure  1.  Plane  Wave  Shielding 


i 


5 


Paramet 


7 


Z.  -  Z  I,  vs  A,„  -  Z  C  Input 

.  in  ws  I  1  ws  T  rr. 

P.  =  -- - ,  v, —  =  -  -----  =  - - - 7T  =  Reflection 


in  Z.  +  Z  V. 

111  WS  -  +  Z 

I  Vs 


+  ^ws  Coefficient 


where  Z^  is  the  input  impedance  to  the  2-port  and  other  circuit  quantities  are 
defined  in  Figure  3. 

Substituting  from  equations  (1)  and  (3)  and  (2)  and  solving  for 


and  — —  yields: 

l2 


=  ~{[l  +  V"S]  cosh0  +  +  --"IslnhO  } 

wL  wL 

=  Reciprocal  of  the  electric  shielding  ratio 

1  Z  z 

-  =  J  {[-^  +  l]  coshO  +  [— —  +  ]  sinhM 

l  2  Zws  ws 

=  Reciprocal  of  the  magnetic  shielding  ratio. 


Then  , 

Magnetic  Shielding  Effectiveness  =  -20  log 


!0  1  l 1 
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Electric  Shielding  Effectiveness  -  -20  log 
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Notice  that  magnetic  and  electric  shielding  ratios  are  identical  if 

Z  -  Z  ,  i  .e.,  for  Z  =  Z  -  Z  ,  the  inverse  ratios  are 
ws  wL  wL  ws  o 

v: =  t:  =  coshG  +  i  +  ?r)slnhe-  (8) 

2  2  o 

By  proper  choice  of  Z  and  Z^g  to  match  the  incident  wave  characteris¬ 
tics,  it  is  possible  to  use  equations  (4),  (5),  and  (8)  as  very  good  approxi¬ 
mations  to  a  number  of  other  more  complex  shielding  configurations.  Physically 
this  is  true  because  the  wavelength  is  typically  much  smaller  within  the  shield 
material  than  outside.  In  most  shield  configurations  shield  thickness  is  small 
compared  with  shield  radii  of  curvature,  and  to  the  characteristics  of  sources 
utilized  in  certain  measurement  structures. 

Multilayer  shields  are  easily  handled  by  using  the  cascading  property  of 
ABCD  parameters  as  shown  in  Figure  4.  Having  obtained  the  effective  ABCD 
parameters  of  the  multilayer  shield,  equations  (4),  (5),  and  (8)  are  used  as 
before  to  evaluate  the  various  shielding  requirements. 


.  -5*  . 
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The  plane  wave  shielding  given  above  is  used  as  a  "reference"  through¬ 
out  this  report.  Each  new  configuration  is,  where  possible,  mathematically 
related  to  the  reference  plane-wave  configuration.  For  example,  as  will  be 
shown,  it  is  possible  to  relate  both  transfer  impedance  and  two-loop/ flat  plate 
data  to  the  above  model. 

The  transmission-line  model  of  isotropic  shields  can  be  extended  to 
handle  plane  waves  incident  at  oblique  angles  [5],  The  required  modification 
is  straight  forward,  and  only  the  characteristic  impedance  and  propagation 
factor  must  be  changed. 

3. 2  Coupling  Mechanisms 

Four  basic  coupling  mechanisms  [6]  degrade  the  performance  of  a  shield. 
The  mechanisms  are  skin  diffusion,  aperture  coupling,  joint  coupling  and  TEM 
penetration . 

Consider  a  solid  (no  apertures,  seams,  etc.)  shield.  An  impinging 
electromagnetic  field  establishes  a  current  density  on  the  exterior  surface  of 
the  shield.  That  current  density  penetrates  by  diffusion  to  the  interior 
surface  of  the  shield  where  a  tangential  electric  field  is  establsiehd.  The 
proportionality  factor  between  the  interior  tangential  electric  field  and  the 
exterior  surface  current  is  called  the  surface  transfer  impedance.  The  overall 
electromagnetic  energy  transfer  by  this  mechanism  is  called  skin  diffusion. 

Skin  diffusion  is  controlled  by  proper  choice  of  material  thickness,  conduc¬ 
tivity  and  permeability. 

Aperture  coupling  occurs  as  the  result  of  an  external  field  exciting 
holes,  cracks,  windows  or  other  openings  in  the  shield.  The  field  established 
in  the  aperture  by  the  external  field  then  couples  to  the  internal  region.  The 
effects  are  typically  difficult  to  calculate  exactly  for  a  variety  of  reasons 
including  strange  aperture  shapes,  complicated  and  changing  interior  and  exte¬ 
rior  environments.  Computer  programs  such  as  IEMCAP  use  shape  approximations 
and  other  idealizations  to  predict  aperture  response.  Aperture  coupling  is 
controlled  by  proper  use  of  screens,  metallic  coatings,  rf  gaskets,  etc. 

Well  formed  joints  in  shield  materials  (i.e.,  joints  of  uniform  construc¬ 
tion  and  good  electrical  contact  without  large  apertures  or  cracks)  may  still 
provide  a  coupling  mechanism  via  the  change  in  surface  impedance  experienced  in 
crossing  the  joint.  Such  joints  usually  are  described  in  terms  of  a  distributed 
joint  transfer  admittance  per  unit  of  joint  width. 

TEM  penetration  is  the  result  of  induced  currents  propagating  along 
wires,  cables,  antenna  feeds  or  other  conducting  paths  into  the  shielded  area. 
These  currents  are  transmission-line-like  currents  (hence  the  name  TEM).  This 
mode  of  coupling  can  cause  further  problems  by  the  induced  currents  generating 
a  field  along  the  wire  which  may  at  a  distant  point  induce  other  troublesome 
currents  in  some  circuit  thought  to  be  well  shielded. 

3. 3  Two-loop/ Inf ini te  Flat  Plate  Configuration 

This  configuration  (Figure  5),  with  the  two  loops  parallel  to  the  flat 
plate  having  a  common  axis  and  spaced  such  that  the  flat  plate  lies  in  the  near 
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field  of  the  loops,  is  closely  related  to  several  popular  and  useful  configur¬ 
ations  for  measuring  the  shielding  effectiveness  to  low- impedance  waves.  The 
low- impedance  (magnetic  field  dominates)  waves  as  viewed  at  the  flat  plate 
result  from  being  in  the  near  field  of  the  loops.  Note  that  the  wave  impedance 
of  the  impinging  wave  varies  over  the  surface  of  the  plate  as  a  function  of 
distance  from  the  source.  If  the  coil- to-plate  spacing  is  a  small  fraction  of 
a  wave  length,  most  of  the  energy  is  concentrated  within  a  narrow  zone  under 
the  coil  so  that  in  approximate  representations  a  constant  source  wave  impe¬ 
dance  may  be  used  with  little  error. 


Shielding  data  from  configurations  of  this  type  for  fiber-reinforced  . 
laminates  have  simple  interpretations  only  for  materials  which  are  essentially 
isotropic  in  the  plane  of  the  flat  plate,  e.g.,  for  multi-layer,  mixed  orien¬ 
tation  graphite  laminates  but  not  for  unidirectional  samples.  This  restriction 
is  a  consequence  of  the  need  for  roughly  circular  currents  to  flow  in  the  plane 
of  the  flat  plate  if  the  incident  magnetic  field  (which  is  largely  perpendicu¬ 
lar  to  the  plate)  is  to  be  terminated. 


Configurations  for  which  the  two-loop/infinite  flat  plate  configuration 
serves  as  a  useful  approximate  model  include  a  variety  of  box-like  structures 
[6]  with  the  flat  plate  shield  forming  a  partition  and  undesired  leakage 
between  loops  inhibited  by  the  enclosure.  Such  configurations  are  typically 
used  at  frequencies  below  100  MHz.  Care  must  be  taken,  however,  since  box 
resonances  can  obscure  the  shield’s  properties. 

•  t 

Moser  [7]  and  Bannister  [8,  9]  have  provided  an  integral  equation  solu¬ 
tion  for  the  shielding  effectiveness  of  the  two-loop/infinite  flat  plate 
geometry  assuming  uniform  current  in  the  loops.  It  is  assumed  that  the  shield 
is  a  good  enough  conductor  that  displacement  currents  in  the  shield  can  be 
neglected.  The  complete  expression  for  shielding  effectiveness  as  given  by 
Bannister  is  given  in  Appendix  A. 


As  shown  in  Appendix  A,  if  it  is  assumed  that  r 1  < 


lair 


T  r ' 
r 


10, 


t  r  * 
r 


20  ’  Trt  >  2’ 


U. 


>  10,  z  »  t  and  z  a,  where  the  quantities  are  defined  in 


the  appendix,  then 


S 


E. 


dB 


8. 686/2  T  t  -1  20  log-^ 
r  10 


T  z 

r 


8.485 


(9) 


This  equation  can  be  shown  to  be  of  the  same  form  as  the  plane  wave 
shielding  equation  provided  an  appropriate  near-field  value  is  used  for  the 
source  and  load  wave  impedances.  Equation  (8)  for  plane  wave  shielding  can  be 
cast  in  the  well-known  Schelkunoff  form  [  2  ]  using  the  following  parameter 
def ini tions 


k 


k  -  1 
k  +  1 


P 


4k 

u  +  l)2 


Using  the  above  parameters  equation  (8)  can  be  written  as; 


(10) 
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Inverse  Shielding  1  -20.  0 

D  =  -  (1  -  qe  )e 

Ratio  p 


(11) 


Shielding  effectiveness  is  then: 


S.E 


dB 


20  lo810 


(12) 


For  the  present  situation,  using  =  jcop^r^  anc*  recognizing  |k|  >>  1,  it 
follows  that  q  -  1  and  —  -  k/4.  Furthermore,  for  T  t  >  2,  the  term 

I  -201  P  r 

20  log  | 1  -  qe  I  is  negligible.  Thus,  equation  (12)  becomes 


S.E 


dB 


-  8.686^  X  t  +  20  log 

r  10 


r  T  r_  >, 
— r-1-. 
2.828y 

^  r; 


(13) 


Following  Moser,  to  compensate  for  the  fact  that  the  near  field  charac¬ 
teristic  wave  impedance  is  actually  not  constant  over  the  shiled,  let  r^  =  z/3. 
Then  equation  (13)  becomes  identical  with  the  simplified  Moser  formula  of 
equation  (9) .  This  indicates  that  under  a  class  of  important  measurement 
conditions  the  plane-wave  shielding  equations  with  appropriate  source  and  load 
wave  impedances  yield  excellent  results.  The  transmission  line  analogy  thus 
applies  to  the  configurations  discussed  in  Sections  1.3.1  and  1,3.2.  Frequent¬ 
ly  a  symmetric  arrangement  with  =  r ^  =  z/2  is  utilized  in  measurements. 

3.4  Quasistatic  Shielding  Formulas  for  Electrically  Thin-Shell  Ellipsoids 


The  boundary  value  problems  for  certain  ellipsoidal-shell  shields 
geometries  as  shown  in  Figure  6  have  been  solved  and  the  corresponding  magnetic 
shielding  effectiveness  calculated  [9].  Following  King,  formulas  for  each  of 
the  ellipsoidal  and  degenerate  ellipsoidal  shielding  formulas  can  be  obtained 
in  the  same  form  as  the  plane  wave  shielding  equations.  These  formulas  are 
useful  in  interpreting  measured  data  from  flat-plate  and  quadrax  structures  as 
well  as  spheres  and  closed  cylinders.  Using  the  notation  developed  for  the 
plane  wave  shielding  to  low- impedance  impinging  waves  can  be  placed  in  the 
form 


Inverse  Magnetic 
Shielding  Ratio 


M 

cosh  0  +  —  sinh  8 

2n 


(14) 


A  similar  relationship  can  be  derived  for  high  wave-impedance  impinging 
signals.  The  equation  is  of  the  same  form  as  the  above  case  but  with  a 
different  wave  impedance  for  the  impinging  signal. 


Inverse  Electric 
Shielding  Ratio 


cosh  8  +  —  sinh 

2n 


(13) 


Equation  (15)  is  different  from  that  given  by  Boeing  [  6]  but  reduces  to  the 
Boeing  form  as  a  special  case. 
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Figure  6.  Quasistatic  Ellipsoidal  Shields. 
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3.5 


Surface  Transfer  Impedance  and  Effective  Conductivity 


Surface  transfer  impedance  has  been  used  for  many  years  as  a  measure  of 
shielding  effectiveness.  Combined  with  other  "two-port"  parameters  it  can  also 
be  used  to  characterize  shielding  materials  in  computer-aided  analysis  programs 
for  determining  interior  and  scattered  fields  of  complex  geometrical  structures 
[10].  For  shields  which  are  thin  compared  to  the  radii  of  curvature  of  the 
shield  and  for  which  wavelength  within  the  shield  is  much  smaller  than  that 
external  to  the  shield,  the  electromagnetic  behavior  of  the  shield  is  essen¬ 
tially  a  local  phenomena.  Each  local  region  may  then  be  considered  planar  [2], 
For  a  planar  shield,  the  two-port  parameters  are  given  in  Figure  7. 

The  impedances  Z^  and  ^21  are  known  as  sur^ace  transfer  impedances 
since  they  relate  field  values  at  opposite  interfaces  of  the  shield.  Surface 
transfer  impedance  is  frequently  measured  using  triaxial  or  quadraxial  config¬ 
urations  [6]  and  the  data  are  reduced  using  the  Schelkunoff  theory  [2].  Notice 
that  for  electrically  thin  samples  (i.e.,  small  0)  cosh  0-1  and 
Zu  -  Z 22  ~  Z  =  ^21  =  ^  csc^  Thus,  for  electrically  thin  shields  a 

measurement  of  the  transfer  impedance  completely  describes  the  shield.  For 
electrically  thick  samples  and  Z  are  not  approximately  equal  to  ar*d 

additional  measurements  are  required  to  totally  characterize  the  shield. 

From  Figure  7  surface  transfer  impedance  written  with  the  new  symbol  Z 
becomes : 


Z  =  ri  csch  0  (16) 

Surface  transfer  impedance  can  be  related  to  the  two-loop/flat-plate 
configuration  through  use  of  the  approximations  for  Equation  (9)  or  the 
equivalent  Equation  (14)  as  follows 


Inverse  Magnetic 
Shielding  Ratio 


cosh  0  + 


2n 


si nh  G 


Substituting 


sinh  9 

n 


1 

rj  csch  0 


“ —  yields 
tr 


R  -  cosh  0  + 


Thus , 


Z 

tr 


Z 


M 


2(R  -  cosh  0) 


(17) 


(18) 


(19) 


For  good  shields  cosh  0  -  1  at  low  frequencies  where  0  is  small  (i.e.,  the 
shield  is  electrically  thin)  and  at  high  frequencies  R  is  much  larger  than 
cosh  0.  Thus,  for  reasonably  good  shields  Z  can  be  written  in  terms  of  the 
shielding  effectiveness  measured  in  the  two-loop/ flat-plate  configuration  as 
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where  z 


=  n  coth  0 


11  *  z22 

Z12  *  Z21  =  n  csch  9 

e  =  yt 

r)  »  wave  impedance  of  shield 
Y  -  propagation  factor  of  shield 

Two-port  Impedance  Representation 

-V 

Note:  H  and  E  directions  assumed  cause 

propagation  toward  the  shield  from 
both  left  and  right  side  of  the 
shield 


l 


i 


(20) 


Ztr  ~  2 (R  -  1) 

where  R  is  obtained  from  the  shielding  measurement  as 

tea 

R.  10 '  20 

At  frequencies  where  the  shield  is  electrically  thin  so  that  ft  is 
sufficiently  small  that  sinh  0  can  be  replaced  by  0,  it  is  possible  to  deter¬ 
mine  an  effective  conductivity  directly  from  Z  in  a  very  simple  manner. 

Such  scalar  conductivity  numbers,  as  mentioned  earlier,  seem  meaningful  for 
multi-layer,  mixed  orientation,  graphite  laminates,  but  not  for  unidirec¬ 
tional  graphite  laminates.  Assuming  a  scalar  effective  conductivity  a  in 
the  plane  of  the  laminate  and  recognieing  that  for  reasonably  good  conductors 

Y  ~  (1  +  j)  v%fuaeff  =  aeffn, 

then 


The  form  of  the  approximation  given  in  Equation  (23)  is  in  excellent 
agreement  with  measured  data  given  by  Boeing  [6]  with  varying  thickness. 
Calculated  conductivities  seem  to  be  in  agreement  with  results  from 
measurements  by  other  techniques  discussed  in  Sections  2.0  and  3.0. 
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As  an  example  of  the  correspondence  between 
ments  and  surface  transfer  Impedance  measurements, 
5208  graphite  samples  measured  by  Boeing  [6],  The 
layups  (0°/45°/90°) . 


two-loop/ f lat-plate  measure- 
consider  the  24-ply  T-300/ 
samples  were  cross-ply 


At  1  MHz  the  measured  magnetic  shielding  from  the 


configuration  is  M.S.E.  =  16  dB. 
plate  spacing  of  1  inch  and  |z^| 


Thus,  R 
=  |  jcoy^b  | 


10 


0.8 


6.31 


0.201  ohm. 


two-loop/ flat-plate 
and  for  a  loop-to- 


Then 


rL 


M 


tr  1 


_  =  0.201 

2 (R  -  1)  2(5.31) 


1.9  x  10  ^  ohm 


as  calculated  from  the  two-loop/flat-plate  measurement  data.  The  corresponding 
1  MHz  direct  measured  value  of  surface  transfer  impedance  is 

|z  I  ~  1.8  x  10  ^  ohm 
tr 


as  measured  in  the  quadrax  configuration.  The  agreement  is  excellent. 
For  the  same  material  at  1  MHz 


11  4 

a  ff  •  yr — r~  =  — - y -  - — - - y-  -  1.73  x  10  mhos/m 

e  1  tr1  (1.8  x  10  )(24  x  5.25  x  10  x  2.54  x  10  ) 

for  T-300/5208  Cross-ply  layup  at  1  MHz. 

Similar  calculations  for  12-ply  HTS/5208  graphite  show  excellent  agree¬ 
ment  between  flat-plate  and  quadrax  data  and  indicate 

O  ---1.5x10  mhos/m 
ef  f 

for  HTS/5208  cross-ply  layup  at  1  MHz. 

3.6  Transverse  Flat  Plate  Samples  in  Waveguide  and  Transmission  Line 
Structures 


Transmission  loss  and  phase  measurements  on  a  flat  plate  sample  com¬ 
pletely  filling  the  transverse  section  of  a  waveguide  or  transmission  line 
structure  can  be  utilized  to  characterize  a  material  elect romagnetically  [11]. 
Provided  the  reflections  from  the  sample  are  not  too  large,  conductivity  and 
permittivity  can  be  determined  analytically  from  the  measured  insertion  loss 
and  phase.  Reference  [11]  provides  a  complete  discussion  on  limitations  of 
the  solution  procedure.  Typical  arrangements  for  rectangular  waveguide  and 
coaxial  line  structures  are  shown  in  Figure  8.  The  analysis  of  these  and  other 
structures  can  be  carried  out  simultaneously  through  the  use  of  the  generalized 
transmission  line  analogy  [2]. 
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3.6*1  General  Case 


Following  Schelkunof f *s  procedure  [2],  let  each  section  of  the  measure¬ 
ment  structure  be  represented  by  a  section  of  generalized  transmission  line. 
The  geometry  being  discussed  is  shown  in  Figure  9.  The  characteristic  imped¬ 
ances  of  the  equivalent  transmission  lines  are  interpreted  as  the  wave 
impedance  of  the  actual  structure.  Equivalent  transmission  line  propagation 
factors  are  equal  to  the  corresponding  factor  of  the  real  structure. 

From  Figure  8, 


X  B 

p 

Q 

C  D 

-i„ 

p 

Q 

(25) 


where 


TQ  =  _VQ/203 

e 

=  Yt 

A  =  cosh  0 

B 

■  ZQ2  sinh 

g  sinh  0 

D 

=  cosh  0  . 

Z02 

From  equation 

(25), 

it  can  be 

^  =  ~ZQ2 

Ip  Z^cosh  0  +  Z^sinh  0 


(26) 


Zo.  _ _ 

Vp  Z^cosh  0  +  Z^^sinh  0 


J03 


(27) 


Let  the  total  impedance  at  P  looking  toward  Q  be  equal  to  Zp  (note,  Zp 
includes  effects  of  both  interfaces).  Then, 


2Zp 

Vp  =  tvv.  =  ^  V. 

2Z„ 


01 


I. 


Ip  '  TlIi  '  'z01  +  ZP  ‘‘ 


(28) 


(29) 


where  V^,  1^  are  the  incident  voltage  and  current,  respectively,  and 
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Material  Under  Test 
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Transverse  Disk  of 
Ma t  e rial  Under  T  e  s  t 


(b) 

Figure  8.  Two  Types  of  Measurement  Strucutres 
Included  in  the  Generalized  Analysis 
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Figure  9. 


Generalized  Transmission  Line  Model  ft 
Measurement  Structure.  Mote: 


vp  = 


V.  +  V 
l  r 


xp  = 


I 


20 


z 


p 


02 


Z^cosh  0  +  Z^sinh 
Z02COS^  ®  +  Z^sinh 


0 

0 


(30) 


Replacing  cosh  0  and  sinh  0  by  their  exponential  equivalents  and  sub¬ 
stituting  from  Equations  (26)  to  (30)  into  Equation  (25),  the  following  result 
is  obtained  (after  considerable  algebra). 


T  ^  -  (1  +  Pv,P)(1  ~  Pv,Q)C 

TI  "  I.  “  .  n  ,-20 

1  1  "  Pv,PPv,QG 


(31) 


where  =  overall  current  transmission 


Z01  z02 


V,P  Z01  +  Z02 


"  =  interface  voltage 


Z03  Z02 

■=  interface  voltage 


’•O  Z03  +  Z02 


It  can  be  shown  that  the  overall  voltage 
section  PQ  is 


coefficient  across  the  section  PQ 


reflection  coefficient  at  P 


reflection  coefficient  at  0. 


transmission  coefficient  across 


the 


T 


V 


If  structures  1  and  3  are  identical  then 


(32) 


T  = 


T 


v 


(i 

i 


2  -20 

Pv" 


(33) 


where 


P 


v 


-  Z 

i 


02 

02 


Insertion  loss  is  thus 


I.L.  =  -20  log1()|T 


(34) 


and  the  insertion  phase  is 


Ac(>  =  U  ~  Angle  of  Transmission  Coefficient 


Given  measured  values  of  insertion  loss  and  insertion  phase  for  a 
specified  structure.  Equation  (33)  can  be  solved  for  the  conductivity  and  the 
real  part  of  the  permittivity  of  the  material  under  test.  This  procedure  is 
demonstrated  in  Section  3.6.2  for  a  rectangular  waveguide  measurement  system. 

3.6.2  Rectangular  Waveguide 

As  a  particular  case  of  the  preceeding  analysis,  let  the  transmission 
structure  be  TE..^  rectangular  waveguide  with  region  2  consisting  of  a  section 
of  guide  of  width  "a,"  height  "b"  and  length  nt,M  totally  filled  with  the 
material  under  test.  Further  assume  that  the  material  can  be  represented  as  a 
lossy  dielectric  with  p  =  p  ,  c  =  e  e’  (1  -  jtan6).  Tan6  is  the  loss  tangent 
of  the  material  and  c1  is  tfie  real  part  of  the  permittivity.  Regions  1  and  3 
are  air-filled  and  have  the  same  width  and  height  as  section  2. 


"TE 


v 


A  2 
!-(-/) 

c 


=  Wave  impedance  of  regions  1  and  3 


z2  = 
TE 


Kr  a 

r  c 


where  n 


•Vf 

■  r.;  (1  -  j  tanfi),  -  -J 


-  -  Wave  impedance 
2  of  region  2 


(36) 


(37) 


-  wavelength  in  freespace  at  measurement  frequency 
A  =  cutoff  wavelength  (e.g.  2a  in  air-filled  waveguide) 


Also, 


z}_  -  z2 

p  =  TE - TE_ 

v  z^;  +  z2 

TE  TE 


(38) 


G  =  yt  =  jw/  p  cVT  -.Aq'I  . 

o  o  r  (TH)  t, 


(39) 


Substitution  into  Equation  (33)  yields  the  desired  result. 
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Given  measured  values  of  insertion  loss  and  phase,  Equation  (33)  can 
now  be  solved  for  and  tan  6.  The  conductivity,  a.  is  related  to  the  loss 
tangent,  tan  6  ,  by  a  =  tan  S.  An  iterative  numerical  procedure  [12]  is 

used  to  solve  Equation  (33).  For  relatively  high  conductivity,  materials  such 
as  graphite  epoxy,  insertion  loss  and  phase  are  quite  insensitive  to  variations 
in  e* f  as  shown  in  Figure  10.  It  is  then  exceptionally  difficult  to  determine 
accurate  values  of  from  measured  insertion  loss  and  phase  by  solving  Equa¬ 
tion  (33).  If  e*  is  independently  known,  accurate  conductivity  values  can  be 
obtained. 

4.0  Scattering  Parameter  Techniques 

The  ease  with  which  scattering  parameters  can  be  measured  [4,  13]  makes 
them  especially  well  suited  for  describing  distributed  circuits  and  most 
electromagnetic  propagation  problems.  The  simple  direct  relationship  between 
scattering  parameters  and  transmitted  and  reflected  power  is  an  added  advantage 
in  interpretiftg  results. 

4.1  Definitions 

Generalized  scattering  parameters  have  been  defined  by  K.  Kurokawa  [2, 
14].  These  parameters  describe  the  interrelationships  of  a  new  set  of 
variables  (a^,  b^) .  The  variables  a.  and  b.  are  normalized  complex  voltage 
waves  incident  on  and  reflected  from1the  itn  port  of  the  system.  They  are 
defined  in  terms  of  the  terminal  voltage  V.,  and  an  arbitrary  reference  imped¬ 
ance  Z.,  as  follows:  1 

l 

V.  +  Z.I. 

i  l  l 


V.  -  Z.I. 

b  ±-  — — hr  <4i: 

2/pTzTf 

where  the  asterisk  denotes  complex  conjugate. 

For  most  measurements  and  calculations  it  is  convenient  to  assume  that 
the  reference  impedance  Z.  is  positive  and  real.  For  the  remainder  of  this 
report,  all  variables  and  parameters  in  scattering  analysis  will  be  referenced 
to  a  single  real  impedance  Rq.  Other  normalization  schemes  are  useful  in 
some  cases.  Generalized  normalization  is  discussed  in  the  literature  [2,  14]. 

The  quantities  used  in  defining  S-parameters  for  a  2-port  system  are 
shown  in  Figure  11.  The  independent  variables  a  and  a^  are  normalized  inci¬ 
dent  voltages  (or  the  corresponding  analog  quanitity) ,  as  follows: 


V,  +  I,  R. 


1  1  0  voltage  wave  incident  on  port  1  =  V  //rTT  (42) 

1 - — - -  il  0 
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(42) 


-  ^2^0  _  voltage  wave  incident  on  port  2 

2'  ^  ^ 


=  lii 

Dependent  variables  and  b^  are  normalized  reflected  voltages: 

y  _  X  g  Voltage  wave  reflected  ^ 

^  _  _JL _ 1  0  _  (or  emanating)  from  port  1  _  r  1 

i’  ^ 


V.  -  I„R, 


Voltage  wave  reflected 


V 


^  _  _2 _ 2  0  _  (or  emanating)  from  port  2  _  r2 

2=  ^ 

The  linear  equations  describing  the  2-port  are  of  the  form: 


S  a  or 


bl  Sllal  +  S12a2 


h?.  S21al  +  S22a2 


The  S-parameters  S  ,  S21  S22  are: 


11 


Input  reflection  coefficient  with  the 

=  output  port  terminated  by  a  "matched" 

load  (Z  =  R  sets  a  =  0)  . 
a  1*  0  2 


S  =  — 
22  a2 


Output  reflection  coefficient  with  the 

=  input  port  terminated  by  a  "matched" 

_n  load  (Z  =  R  and  V  =  0). 
a^  U  S  U  S 


21 


a2=0 


Forward  transmission  (insertion)  gain 
=  (or  loss)  with  the  output  port 
terminated  in  a  "matched"  load. 


12 


ar° 


Reverse  transmission  (insertion)  gain 
(or  loss)  with  the  input  port 
terminated  in  a  "matched"  load. 


(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 

(51) 
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The  relationships  between  a^,  a^,  b^,  b^  and  the  various  power  waves  of 
interest  are  given  by 

2 

|a^|  =  Power  incident  on  the  input  port 

-  Power  available  from  a  source  of  impedance 

2 

| |  -  Power  incident  on  the  output  port 

=  Power  reflected  from  the  load 

|b^|  =  Power  reflected  from  the  input  port 

=  Power  available  from  a  source  minus  the  power  delivered  to 
the  input  port 

i  1 2 

b  =  Power  reflected  or  emanating  from  the  output  port 

2 

=  Power  incident  on  the  load 


=  Power  that  would  be  delivered  to  a  R^  load 

Hence  S-parameters  are  directly  related  to  power  gain  and  mismatch  loss, 
quantities  which  are  often  of  more  interest  than  the  corresponding  voltage 
functions : 


*  i 2  _  Power  reflected  from  the  input  port 

*  11*  Power  incident  on  the  input  port 

|s  | 2  Power  reflected  from  the  output  port 

*  22*  Power  incident  cn  the  output  port 

Power  delivered  to  a  Rn  load 

Is  \  =  — — — ■  -  — ■  -  - - - - 

1  21 1  Power  available  from  a  R^  load 

=  Transduce  power  gain  with  R^  load  and  source 

2 

=  Reverse  transducer  power  gain  with  load 
and  source. 

4 . 2  System  Calculations  with  Scattering  Parameters 

Scattering  parameters  are  particularly  convenient  in  calculating  trans¬ 
mitted  and  reflected  power  or  related  quantities  such  as  shielding  ef feet iveness . 
The  transfer  parameters  S  and  ate  a  measure  of  the  complex  insertion  gain, 
and  the  driving-point  parameters  S-  and  are  a  measure  of  the  input  and 
output  mismatch  (i.e.,  reflection)  loss.  As  dimensionless  expressions  of  gain 
and  reflection,  the  parameters  not  only  give  a  clear  and  meaningful  physical 
interpretation  of  the  system  performance  but  also  form  a  natural  set  of  param¬ 
eters  for  use  with  signal  flow  graphs  [15].  It  is  not  necessary  to  use  signal 
flow  graphs  in  order  to  use  S-parameters ,  but  flow  graphs  greatly  simplify 
S-parameter  calculations. 


27 


In  a  signal  flow  graph  each  port  is  represented  by  two  nodes  (one  for 
each  variable) .  Node  a  represents  the  wave  coming  into  the  system  at  port  n 
and  node  b  represents  Phe  wave  leaving  the  system  at  port  n.  The  complex 
scattering  coefficients  are  then  represented  as  multipliers  on  branches 
connecting  the  nodes  within  the  system  and  in  adjacent  systems.  Figure  12 
summarizes  the  above  statements  and  presents  the  flow  graph  representation  of 
the  2-port  system  shown  in  Figure  11. 

The  simplification  of  system  analysis  by  signal  flow  graphs  results 
from  the  application  of  the  "non-touching  loop  rule"  as  described  in  [15]  and 
summarized  in  Figure  13.  A  first  order  loop  is  defined  as  the  product  of  the 
branches  encountered  in  moving  from  a  node  in  the  direction  of  the  arrows  back 
to  that  original  node.  A  second  order  loop  is  defined  as  the  product  of  any 
two  non-touching  first  order  loops.  Nontouching  loops  have  no  nodes  in  common. 
An  n^h  order  loop  is  defined  as  the  product  of  any  n  non-touching  first  order 
loops.  This  rule  applies  the  generalized  formula  of  Figure  13  to  determine  the 
transfer  function  between  any  two  nodes  within  a  system.  The  non-touching  loop 
rule  is  applied  to  calculating  the  transducer  power  gain  of  a  2-port  system  in 
Figure  14. 

Cascading  of  2-port  systems  is  easily  accomplished  through  use  of  signal 
flow  graphs  and  the  S-parameter  definitions  given  in  Equations  (46-51).  The 
cascading  is  easiest  if  done  2  systems  at  a  time.  This  process  and  the  result¬ 
ing  combined  S-parameters  are  illustrated  in  Figure  15. 

4.3  Parameter  Relationships  for  Two-Ports 

It  is  convenient  at  times  to  be  able  to  convert  between  the  various 
parameter  representations  for  2-port  networks.  The  equations  below  provide 
those  conversion  relations  most  useful  in  shielding  analysis.  Other  conversions 
are  in  the  literature  [3,  4], 

4.3.1  Z  to  S 


Given  a  normalized  impedance  matrix 


Z 

Z 


11  Z12 
21  Z22 


(52) 


the  corresponding  S-matrix  is 


4— > 

s 


s 

s 


11  S12 
21  S22 


(53) 


where 
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•  EACH  VARIABLE  BECOMES  A  NODE. 

•  EACH  PARAMETER  BECOMES  A  BRANCH. 

•  BRANCHES  ENTER  DEPENDENT  VARIABLE  NODES  AND  EMANATE  FROM 

INDEPENDENT  VARIABLE  NODES. 

•  EACH  NODE  IS  EQUAL  TO  THE  SUM  OF  THE  BRANCHES  ENTERING  THAT 

NODE. 

b1  *S11  a1 +s12a2  b2  =  S21  a1  +s22a2 


S-J2  a2 


Figure  12.  Flow  Graph  Notation 
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TWO-PORT 

NETWORK 


(Cont ' d) 


NONTOUCHING  LOOP  RULE 

T  =  P,  (1  -  ZL  +  Z  L  (2)(1)  -  Z  L  (3)(1)  +  •  •  •  ]  +  P2  (  1  -  ZL  (1>(2)  +  •  •  • 

1  -ZL(1)  +  ZL  (21  —  Z  L  (3|  +•  •  •  • 


WHERE 

ZL(1)  SUM  OF  ALL  FIRST  ORDER  LOOPS. 

IL(n)  SUM  OF  ALL  nth  ORDER  LOOPS. 

Pi;  P2;  etc.  PATHS  CONNECTING  VARIABLES  IN  QUESTION. 

IL(1)(1)  SUM  OF  THOSE  FIRST  ORDER  LOOPS  WHICH  DO  NOT  TOUCH  P, . 

£L  (n)m  SUM  OF  THOSE  nth  ORDER  LOOPS  WHICH  DO  NOT  TOUCH  Pm. 

T  RATIO  OF  DEPENDENT  TO  INDEPENDENT  VARIABLE. 


INPUT  IMPEDANCE  OF  A  TWO-PORT  WITH  ARBITRARY  LOAD 


1  “  1 


lid  -s22  rL 


p2  =  s21  rLs12 


b1  S11  l1  "  s22  1  L*  +  S21  1  L  S12 

»i  '  i-s22rL 


S21  S12  1  L 

sn  _sn  +  i  -s22rL 


Figure  13.  Nontouching  Loop  Rule  for  Use  in  Evaluating 
Transfer  Functions 
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I  ^3  l 

P-ITsI2] 


b2  _  _ S21 

bs  1_S1 1  !'s  "  s22  rL  ~  S21  S12  r  L  rs  +  S11  rs  S22  rL 

g  =  ls2i  1  2  H  ~  1  rsl2  >  <1  -  I rL I 2  > 

°r  Id  -  sn  rs »  (i  -  s22  rL  )  -  s21  s12  rL  rs|2 

Figure  14«  Transduce  Power  Gain  Calculations 


32 


4.4 


Multiport  Sv s_t  oms 


Scattering  parameter  analysis  is  easily  extended  to  multiport  systems. 
For  a  4-port  system  the  scattering  equations  would  be: 


where 


S 

S 

s 

s 


11 

S12 

S13 

S14 

21 

S22 

S23 

S24 

31 

S32 

S33 

S34 

41 

S42 

S43 

S44 

(56 ) 
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a=a=a  =0 
2  3  4 


Input  reflection  coefficient 
at  port  1  with  matched  load  at; 
all  other  ports. 


21 


a2=a3=a4=0 


Forward  transmission  (insertion) 
gain  (or  loss)  with  ports  2,  3, 

4  all  terminated  in  matched  loads. 


etc . 

The  extension  to  n-ports  is  straightforward.  Multiport  techniques  are  necessary 
in  treating  shielding  effectiveness  of  anisotropic  shields. 

5 . 0  An  Anisotropic  Model  for  Plane  Wave  Propagation  (with  Applications _ t o 

Single  and  Multilayer  Advanced  Composites) 

Advanced  composite  materials  are  laminates  made  up  of  a  number  of  indi¬ 
vidual  layers  bonded  together.  Each  layer  consists  of  a  unidirectional  array 
of  long  fibers  embedded  in,  and  firmly  bonded  to  a  matrix.  The  basic  building 
blocks  of  any  specific  composite  are  the  types  of  fibers  and  matrix  involved. 
Some  fiber/matrix  systems  are:  boron/epoxy,  graphite/epoxy,  Kevlar /epoxy , 
graphi te/polyimide  and  graphi te/thermoplastic .  Other  materials  contain  grids 
of  various  conductive  wires. 

The  purpose  of  this  section  is  to  describe  a  macroscopic  electromagnetic 
model  for  plane  wave  propagation  through  fiber-reinforced  and  related  materials. 
Ft  is  convenient  to  distinguish  between  "ply"  and  "layer.”  A  ply  is  the  basic 
composite  material  unit.  A  layer  is  made  up  of  plies  with  uniform  fiber 
orientation.  Thus,  both  plies  and  layers  are  unidirectional  but  a  layer  may 
consist  of  more  than  one  ply.  In  a  multilayer  structure,  fibers  in  adjacent 
layers  have  different  orientations.  A  multilayer  flat  panel  with  arbitrarily 
polarized  incident  plane  waves  is  considered.  The  fibers  in  each  layer  may  be 
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oriented  in  any  desired  direction  permitting  both  unidirectional  and  mixed- 
orientation  samples  to  be  considered. 

On  a  microscopic  basis  advanced  composites  are  clearly  both  inhomogene¬ 
ous  and  anisotropic.  From  a  macroscopic  point  of  view,  the  fiber  spacings  are 
sufficiently  close  that  over  a  wide  frequency  range  (perhaps  dc  -  18  GHz)  it 
appears  that  the  inhomogeneities  can  be  averaged  out,  at  least  to  first  order. 
The  anisotropic  nature,  however,  of  a  unidirectional  sample  must  be  taken  into 
account.  For  mixed  orientation,  multilayer  samples  the  anisotropic  effects  are 
important  unless  the  layers  are  so  thin  electrically  that  considerable  layer- 
to-laver  averaging  occurs. 

Based  on  the  above  discussion,  it  appears  reasonable  to  model  each  layer 
as  a  Mquasi-homogeneousM  anisotropic  material.  Both  fiber  and  matrix  materials 
currently  under  consideration  are  nonmagnetic  so  that  the  permeability  of  these 
composite  materials  can  be  taken  to  be  essentially  p  .  Both  permittivity  and 
conductivity  parameters  may  be  anisotropic.  It  is  assumed  that  the  principal 
axes  of  the  permittivity  and  conductivity  properties  are  the  same  so  that  both 
the  permittivity  and  conductivity  tensors  may  be  diagonalized  simultaneously. 

As  a  labor  saving  device  the  usual  combination  of  conductivity  and  permittivity 
into  a  single  complex  permittivity  tensor  is  assumed.  Two  sets  of  coordinates 
are  utilized  as  shown  in  Figure  16.  The  x,  y,  z  coordinates  are  the  "propaga¬ 
tion  coordinates,"  and  all  the  final  equations  are  expressed  in  terms  of  xyz- 
components.  The  ^  -3  coordinates  are  called  "material  coordinates."  The 

material  coordinate  axes  are  aligned  with  the  principal  directions  of  the 
composite  material. 


Analytically  it  is  possible  that  the  "material  a^es"  ,  £  within 

each  layer  be  skewed  with  respect  to  each  of  the  "propagation  coordinate"  axes 
x,  y,  z.  Physically,  however,  a  laminate  is  made  up  of  parallel  layers  with 
the  fibers  of  adjacent  layers  lying  in  parallel  planes  but  having  different 
fiber  orientation  wi th in  the  plane,  i.e.,  the  material  coordinate  ^  will 
usually  coincide  with  z  in  every  layer  of  a  multilayer  sample,  while  the  direc¬ 
tion  of  £  ,  ^2  with  respect  to  x,  y  varies  from  layer  to  layer. 

5.1  Unidirectional  Samples 


Consider  a  unidirectional  flat  plate  of  advanced  composite  material. 

The  principal  directions  (along  the  fibers,  perpendicular  to  the  fibers  but 
parallel  to  the  plane  of  the  layer,  and  perpendicular  to  both  the  fiber  direc¬ 
tion  and  the  plane  of  the  layer)  are  used  to  define  "material"  coordinates. 

Let  be  parallel  to  the  fibers,  C ^  perpendicular  to  the  fibers  but  parallel 
to  the  plane,  perpendicular  to  both  fibers  and  the  plane  of  the  layer.  In 
these  coordinates  the  complex  permittivity  matrix  is  diagonal  with  each  of  the 
three  components,  in  general,  different. 
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Figure  16,  Propagation  (x,  y,  z)  and  Material  ^2’  C°or<3inates  ^or 

Treating  Anisotropic  Material  with  ’’OpticAxis"  in  z  Direction. 
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/V  O.  . 

where  e..  =  c.  -  i  -  . 

11  11  (jl) 

■4 — ► 

In  terms  of  ^3  the  complex  permittivity  tensor  £  is  diagonal. 

a  i  i 

The  individual  elements  of  £,  namely  £..  =  £.,  -  i  ,  can  be  measured  or 

J  ii  ii  J  a)  * 

calculated  for  a  single  unidirectional  sample  of  each  fiber/matrix  combination. 
The  matrix  in  each  case  is  normally  a  good  dielectric  while  the  fibers  vary 
from  modest  conductors  (graphite)  to  poor  dielectrics  (boron)  to  good  dielec¬ 
trics  (Kevlar).  Thus,  the  relative  importance  of  £..  and  a.,  in  the  complex 
permittivity  £  is  a  function  of  the  fiber/matrix  combination  and  frequency. 

5.2  Wire  Grid  Model 


Prior  to  detailed  analysis  of  the  rather  complicated  anisotropic  model 
for  fiber-reinforced  materials,  it  is  instructive  to  examine  some  less  compli¬ 
cated  geometries  [16]  to  perhaps  determine  some  of  the  dominant  features  of 
shielding  to  be  expected  from  such  materials.  Consider  an  array  of  infinitely 
long,  identical,  parallel,  perfectly  conducting  wires  as  shown  in  Figure  17. 
This  is,  admittedly,  a  very  crude  model,  but  it  is  useful  in  illustrating  an 
important  point.  Let  the  wires  have  a  diameter  D,  and  a  spacing,  S,  and  let  a 
plane  wave  having  wavelength.  A,  be  incident  normally  on  the  grid.  The  inci¬ 
dent  electric  field  may  be  polarized  either  parallel  (i.e.,  E  )  or  perpendicu¬ 
lar  (i.e.,  E  )  to  the  axis  of  the  wires.  The  equivalent  circuit  of  the  grid 
as  seen  by  the  incident  wave  is  given  by  Marcavitz  [17]  and  is  shown  in  Figure 
18.  When  S/A  <<  1  the  circuit  parameters  are: 


where 
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Incident 

Field 


(a)  Front  View 


(b)  Top  View 


Figure  17.  Wire  Grid  Approximation  of  a  Unidirectional 
Fiber-Reinforced  Material  (After  Bodnar 
[16]) 


(a)  For  Parallel 
Polarization 


(b)  For  Perpendicular 
Polarization 


Figure  18.  Equivalent  Circuit  of  Wire  Grid  at  Normal 
Incidence.  (After  Bodnar  [16]) 
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and  Zg  and  V  are  the  characteristic  impedance  and  admittance  of  free  space, 
respectively. 

Inspection  of  (58)  and  (59)  shows  that  X  and  approach  zero  as  A 
becomes  large,  indicating  that  the  inductor  of  figure  l8a  shorts  the  trans¬ 
mission  line  at  low  frequencies  and  little  power  is  transferred  through  the 
wire  grii.  Thus,  the  grid  acts  as  an  effective  shield  to  parallel  polariza¬ 
tion  at  low  frequencies.  Inspection  of  (60)  and  (61)  shows  that  B  approaches 
zero  and  ‘hat  B^  becomes  large  as  A  becomes  large.  This  indicates  that  the 
shunt  capacitors  in  Figure  18b  act  as  open  circuits  and  the  series  capacitor 
acts  as  a  short  as  A  becomes  large.  This  corresponds  to  a  large  amount  of 
energy  being  transferred  through  the  grid  and  shows  that  the  grid  is  not  a 
good  shield  for  perpendicular  polarization. 

How  long  do  the  wires  have  to  be  in  order  to  obtain  good  shielding 
characteristics  for  parallel  polarization?  What  would  the  shielding  charac¬ 
teristics  be  if  the  infinitely  long  wires  of  Figure  17  were  replaced  by  a  two 
dimensional  array  of  short  nontouching  wires?  The  shielding  characteristics 
now  turn  out  to  be  bad  for  both  polarizations  as  will  be  illustrated  below. 

The  conclusion  clearly  is  that  a  large  number  of  noncontacting  fibers  will  not 
provide  good  shielding.  Long  conducting  paths  (compared  to  wavelength??)  are 
necessary  tor  good  shielding. 

A  model  for  fibers  consisting  of  short,  noncontacting  plates  is  shown 
in  Figure  19.  Analysis  of  such  structures  has  been  carried  out  by  Chen  [18] 
and  Montgomery  [19].  General  equations  are  provided  by  both  Chen  and 
Montgomery,  but  the  equations  must  be  programmed  for  a  digital  computer  to 
obtain  numerical  results.  Sample  calculations  are  presented  by  both  authors, 
but  no  general  design  information  is  given.  The  general  trends  indicate  low 
attenuation  for  both  polarizations  when  the  plate  dimensions  are  small  compared 
to  wavelength.  A  reasonable  conclusion  appears  to  be  that  noncontacting  fibers 
that  are  electrically  small  (i.e.,  much  smaller  than  a  wavelength  in  their 
major  dimension)  do  not  provide  effective  shielding  at  RF  frequencies. 

5 . 3  Unidirectional  Sample  with  Normally  Incident,  Arbitrarily  Polarized 

Plane  Wave 

Assume  an  infinite  plane  sheet  of  unidirectional  material  of  thickness 
t  oriented  perpendicular  to  the  z-direction.  Fibers  are  oriented  in  the 
^-direction  as  shown  in  Figure  16.  The  ^  and  z  axes  are  aligned.  The 
direction  of  propagation  of  the  incident  plane  wave  is  +z.  The  permeability 
V  is  essentially  equal  to  \i  .  Permittivity  and  conductivity  are  combined  in 
the  complex  permittivity  tensor  given  bv  equations  (57).  In  material 
coordinates  there  are  two  possible  normal  modes  for  propagation  in  each  direc¬ 
tion,  +  and  -z.  These  modes  are  not  coupled  to  one  another  and  correspond  to 
polarization  perpendicular  and  parallel  to  the  fiber  direction,  respectively. 
Using  the  transmission  line  analogy  [2]  the  response  to  the  normal  modes  ran 
be  written  in  material  coordinates  as  follows. 
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(b)  Top  View 


Figure  19.  Model  Consisting  of  an  Infinite  Periodic 
Array  of  Thin  Conducting  Plates  on  a 
Dielectric  Sheet.  (After  Bodnar  [16]) 
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For  E  parallel  to  the  fiber  direction  (£^) 


\<o  * 

“coshO^  H^sinhO^  “ 

"Vs + 

H  (C„) 

— sinhG,  coshO- 

H,  (C,  +  t) 

c2  3  -1 

Ln1  1  1  J 

L  C2  3  J 

where 


°l  =  u,/uo  £11  =  V 


(62) 
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For  E  perpendicular  to  the  fiber  direction  but  parallel  to  the  plane 


of  the  sample  (^ ) 
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-H,  (r  ) 

'>1  3  -J 

ti 

coshG^ 


r^sinhQ  * 


— sinh09  coshO 

Ho  ^  2 


L  -H  <C3  +  d)  J 


where  9  »  u>/y  j?  t  »  k0t 

l  o  22  2 


n2  '  /  ^22 

E  ,  H  ,  E  ,  H  are  components  of  a  single  arbitrarily  polarized 

h  '2  4  Q1 

plane  wave.  It  is  convenient  to  collect  these  terms  into  a  single  4-port 
matrix  equation  relating  tangential  fields  at  input  and  output  interfaces. 
E  and  H  terms  are  slightly  rearranged  in  the  new  matrix.  Note,  there  is 
no  coupling  between  the  and  £2  polarizations. 


VV 

vv 


coshG^ 
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n^sinhG^ 

0 

0 

cosh©2 
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n2sinh©2 
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0 

coshO^ 

0 

0 

— sinh0„ 

n2  2 

0 

coshG2 

Er  (?_  +  t  ) 
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E_  (5,  +  t  ) 
^2  3 

H  (5,  +  t ) 

C2  3 


-H  (5,  +  t  ) 
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(63) 
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To  transform  from  components  expressed  in  material  coordinates  to 
those  in  propagation  coordinates,  it  can  be  seen  from  the  geometry  of  Figure 
16  that 


A  =»  A  cos<{>  +  A  sind> 
Ci  x  y  y  ^ 


A  =» -A  sin<(>  +  A  cosd> 
C2  x  y 


where  A  represents  either  an  electric  or  magnetic  field  component.  Substitu¬ 
ting  from  equation  (64)  into  equation  (63)  yields  the  following  equation: 


coshQ ^ 

0 

I 

sinh0„ 


0 

cosh0, 

4 

0 

sinh0,. 


r)^sinh0^ 

0 

cosh0. 


n2sinh6J 


cosh0. 


where 


=  E  (z)cos(J)  +  E  (z)sin4> 
ll  x  y 

E^2  =  -Ex(z)sinc()  +  Ey(z)c.osd) 

=  -H^(z)sin(j)  +  Hy(z)cosc() 

Hi2  =  -Hx(z)cos^  ~  Hy(z)sin<|> 

E  -  =  E  (z+t)ccs<t)  +  E  (z+t)sin<f> 
ol  x  y 

Eq2  =  -Ex(z+t)sin<()  +  Ey(z+t)cos<|> 
=  -Hx(z+t)sin(()  +  H  y(z+t)cos0 
H^2  =  -H^(z+t)cosd>  -  Hy(z+t)sin<{>. 


Algebraically  manipulating  equation  (65)  leads  to  equation  (66)  given  below. 
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E  (z+t) 
y 

(z+t) 

-Hx(z+t) 


43 


where 


2  .  2 
^  =  D^^=coshO^  cos  <p  +  coshG^  sin  <J 

2  2 
^22  =  D2^=C°sh^l  sin  ^  +  cos^2  cos  ^ 
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A2.t 

II 

CN 

a 

ii 

D2.I 
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COS([) 

(ri,  sinhOj  - 
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C12 
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Other  quantities  have  been  defined  earlier.  Notice  that  x  and  y  polariza¬ 
tions  are  coupled  unless  <p  =  0°  or  90°.  Thus,  except  for  these  two 
special  cases  an  x-polarized  incident  wave  will  produce  both  x  and  y 
polarized  reflected  and  transmitted  waves. 


A  still  more  condensed  equation  results  if  matrix  notation  is  further 
exploited . 


"Et(z)l  Ta  b  I  r£t(z+tr 

Jit(z)j  V  eT  J  [h  (z+t^ 


where 


Et(-)  = 


Transverse 

— ► 

Components  of  E  at 
either  z  or  z+t 


v> 


y-> 


Transverse 
Components  of  H  at 
either  z  or  z+t 


A11  A12| 


,  etc . 


and  A^j ,  B  _ ,  coefficients  were  defined  above.  This  4-port 

structure  is  illustrated  schematically  in  Figure  20. 
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Output 

Region 


Output 


5.4 


Multilayer  Samples  with  Normally  Incident,  Arbitrarily  Polarized  Plane 
Waves 


To  analyze  multilayer  structures,  construct 


B 

D 


for  each  layer  taking 


into  account  the  different  ((>  and  G  for  each  layer.  Obtain  the  overall  ABCD 
matrix  as  the  product  of  the  matrices  for  the  individual  layers.  Thus,  for 
a  90°,  45°,  90°  multilayer  sample 


A  B 

A90°  B90o 

<— >  4 — > 

A45°  B45° 

< — >■  ► 
A90°  B90o 

^ —  >  +  —> 

4r~Y  4— ► 

4 — »-  4 — ► 

C  D 

C90°  D90° 

r  n 

45°  45° 

r  n 

90°  90° 

(68) 


The  resulting  ABDC  matrix  relates  the  "total11  (i.e.,  sum  of  incident  and  re¬ 
flected)  fields  at  the  input  side  of  the  sample  to  the  corresponding  quantities 
on  the  output  side.  Most  often  the  desired  end  result  of  such  analysis  is 
various  ratios  of  traveling  wave  (not  "total")  fields.  To  facilitate  calcula¬ 
tions  such  as  determining  power  reflected,  power  transmitted,  power  converted 
from  one  polarization  to  another,  etc.,  it  is  convenient  after  the  overall 
ABCD  matrix  is  computed  to  convert  to  scattering  parameters  [4].  Scattering 
parameters  provide  direct  relationships  between  the  traveling  wave  quantities. 
Clearly,  scattering  parameters  or  their  transmission  counterpart  could  be  used 
to  describe  each  laver  and  then  combined  to  provide  the  overall  multilayer 
matrix.  However,  the  determination  of  the  overall  ABCD  matrix  and  then  its 
conversion  to  scattering  parameters  is  computationally  much  more  efficient  in 
the  class  of  problems  being  considered. 


5 . 5  Scattering  Matrix  for  a  Multilayer  Anisotropic  Structure  with  Normally 
Incident,  Arbitrarily  Polarized  Plane  Waves 


Either  normalized  or  unnormalized  scattering  parameters  [4]  may  be  used. 
In  many  cases  normalized  parameters  are  a  distinct  advantage.  At  this  point 
there  is  little  advantage  to  be  gained  from  normalization  and  some  danger  of 
confusion  arising  from  the  normalization  process.  Therefore,  unnormalized 
scattering  parameters  are  used. 


Let  the  fields  on  the  input  side  of  the  sample  be  expressed  as  follows: 
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where  n  =  intrinsic  impedance  of  input  medium  and  the  subscript  i  refers  to 
incident  traveling  wave  components  while  r  refers  to  reflected  or  reverse 
traveling  wave  components. 


Similarly,  fields  on  the  output  side  of  the  sample  are  given  in  equa¬ 
tion  (70).  Note,  as  is  conventional  in  scattering  analysis  "incident  waves" 
travel  towards  the  sample  while  "scattered  or  reflected  waves"  travel  away 
from  the  sample.  This  means  that  incident  waves  on  opposite  sides  of  the 
sample  travel  in  opposite  directions. 

Exi(z+t)  +  Exr(z+t) 

E  . (z+t)  +  E  (z+c)  (70) 

yi  yr 

V<2+t> +  V2+t>  -,W2+t>l 

-[Hx.(z+t)  +  Hxr(z+t)I  -7^[Eyl(z+t>  -Eyr(z+r>] 

where  =  intrinsic  impedance  of  output  medium. 

If  there  are  no  sources  on  the  output  side  of  the  sample,  the  incident 
signals  from  that  side  will  be  zero,  i.e.,  E  .  (z  +  t)  *  E  . (z  +  t )  =  0 .  It 
is  assumed  that  this  is  true  in  the  f  ollowin^analysis .  Substituting  equa¬ 
tions  (69)  and  (70)  into  equation  (67)  and  solving  for  the  scattered  or 
reflected  components  in  terms  of  the  nonzero  incident  components  yields  the 
following  relationship.  Note  that  since  E  . (z  .+  t)  =  E  .(z  +  t)  =0,  only  the 
first  two  columns  of  the  scattering  matrix^are  requiredX£or  this  analysis. 


E^(z+t)  = 
E  (z+t)  = 

y 

Hy(z+t)  = 

-H  (z+t)  = 
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Equation  (71)  can  be  used  to  calculate  the  quantities  of  interest  in  electro¬ 
magnetic  shielding  problems.  The  procedure  is  to  first  calculate  the  total 
4-port  ABCD  parameters  for  the  multilayer  sample,  as  indicated  in  Section  5.3. 
The  values  of  A..,  B..,  C..,  D..  obtained  are  used  in  equation  (71)  to  deter¬ 
mine  the  overall  scalier in^  parameters  for  the  multilayer  sample.  Figure  21 
and  22  illustrate  shielding  effectiveness  calculated  for  some  single-layer  and 
multilayer  graphite  epoxy  shields.  The  value  of  £f  is  not  critical  for  the 
values  of  conductivity  utilized.  Indeed  varying  f from  1  to  at  ;>ut  100 
has  no  noticeable  effect  on  the  shielding  effectiveness.  Figure  21  illustrates 
the  effect  on  shielding  effectiveness  of  varying  polarization  of  the  wave 
incident  on  a  single  anisotropic  layer.  Data  is  plotted  for  polarization 
parallel  to  the  fiber  direction  (f  -  0°),  perpendicular  to  the  fiber  direction 
($  =  90°)  and  inclined  at  45°  to  the  fiber  direction  (if  -  45°).  As  would  be 
expected,  much  larger  attenuation  is  experienced  when  the  wave  is  polarized 
parallel  to  the  fibers.  Figure  22  illustrates  shielding  effectiveness  for 
4-layer  (0°,  +  45°,  90°)  and  7-layer  (0°,  +  45°,  90°,  +  45°,  0°)  structures. 

The  incident  wave  in  each  case  is  polarized  parallel  to  the  fiber  direction 
of  the  first  layer. 

5 . 6  Scattering  Matrix  for  a  Multilayer  Anisotropic  Shield  with  Incident 

Plane  Wave  of  Arbitrary  Polarization  and  Angle  of  Incidence 

Consider  a  shield  consisting  of  planar,  homogeneous,  anisotropic  layers 
arranged  as  shown  in  Figure  23.  Each  layer  is  assumed  uniaxial  with  the  optic 
axis  parallel  to  the  xy  plane.  Input  and  output  regions  are  homogeneous  and 
isotropic  with  relative  permittivities  c;"  and  respectively.  A  plane  wave 

of  radian  frequency  is  incident  from  the  region  z  <  0  on  the  shield  with 
angle  of  incidence 
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Shielding  (dB) 


Figure  21.  Shielding  effectiveness  of  a  single  layer  anisotropic 

shield  for  various  fiber  orientations.  A  plane  wave 

polarized  in  the  <t>  =  0°  direction  is  normally  incident 

on  the  innut  side  of  the  shield,  shield  parameters  are 

/  ? 

t  =  .00525  inches,  o  =  2  x  10  mhos/m,  -  2  x  10 
mhos/m,  ~  3 . 
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Frequency  (H^) 


Figire  22.  Shielding  effectiveness  of  4-layer  and  7-layer 
anisotropic  shields .  Fiber  orientations  for  4-layer  shield  are 

(0°,  +  45°,  90°);  for  7-layer  shield  are  (0°,  +  45°,  90°,  4-  45° 
Each  layer  has  t  =  .00525  inches,  o  =  2  x  104  mhos/m, 

O  =  2  x  10^  mhos/m,  e  =  e  „  =  3. 


% 

-  Incident 

Wave 

Figure  23.  Multilayer  Anisotropic  Shield  with  Incident  Plane  Wave  at 
Arbitrary  Angle  of  Incidence.  Shield  Infinite  in  Extent 
in  xy-Plane.  Plane  of  Incidence  is  the  xz-Plane. 


To  evaluate  the  shielding  provided  by  this  structure,  the  problem  may 
be  stated  as  follows.  Given  the  parameters  of  the  anisotropic  layers  of  the 
shield  plus  the  frequency  and  direction  of  propagation  of  the  incident  wave, 
determine  the  reflected  and  transmitted  waves.  Using  the  techniques  of 
Section  5.3  of  this  report,  the  above  problem  reduces  to  that  of  determining 
the  transfer  matrix  of  the  stratified  anisotropic  shield  for  the  specified 
excitation.  The  folloiwng  derivation  of  the  necessary  transfer  matrix  is 
similar  to  the  optic  problem  solved  by  Semenko  and  Mironov  [20]. 

Within  each  anisotropic  layer.  Maxwell’s  curl  equations  become 

->  + 

V  X  E  «  -jwu  H 

-  —  °  <->  -  (72: 

V  X  H  =  (a  +  we)  E 

where  e  and  c  are  the  matrix  representations  of  the  layer  conductivity  and 
permittivity,  respectively,  and  unity  relative  permeability  has  been  assumed. 

Eliminating  H  from  (72)  yields  the  wave  equation  for  this  type  of 
anisotropic  media. 

V2E  =  V(V  •  E)  -  B2<e>E  (73' 

O 

where 

w 

Bo  2 
c 

4 — ►  1  4 — > 

<e>  =  £  -  j  -  0  =  Complex  permittivity 

o 

•4 — >■ 

e  =  Matrix  representation  of  relative  permittivity 
=  Permittivity  of  free  space,  and 

4 — ►  4 — y 

£^  and  0  are  assumed  to  have  the  same  optic  axis. 

Based  on  the  assumption  that  each  layer  is  uniaxial  and  recognizing 
that  the  matrix  representation  of  the  complex  permittivity  is  symmetric  [21], 
<£>  can  be  written  in  the  following  form. 
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where 
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Cj  cos2  G  ♦  sin2  c 

.  2  2 
sm  £  cos  c 


e12  =  ^E2  "  EP  Sin  C  COS  C 


e33  =  e2 


are  the  principal  axis  complex  permittivities 

C  is  the  angle  between  the  x-axis  and  the  optic 
axis. 

The  tangential  field  components  and  their  first  derivatives  are  contin¬ 
uous  across  any  interface  parallel  to  the  xy  plane.  Thus, 


f -0 


(75) 


-r—  -  - i 8  7z  sin  $  =  constant  =  -ik 

3x  o  or  o  x 


where  E  is  the  relative  permittivity  of  the  input  medium. 

Substituting  (74)  and  (75)  into  (73)  yields  for  the  z-component  of  E 

-Jk 


E  * 
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~  D  E 
f  z  x 


(76) 


where 

k  -  3  /z~  sin  <j> 
x  o  o  o 

D  =  d/dz 
z 

2  2 

f  =  6  -  V  . 

o  33  x 

The  tangential  field  components  E  and  E  must  satisfy  the  following  homogen¬ 
eous  matrix  equation.  X  ^ 
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S2e  -  k2  +  D2 
o  22  x  z 


(77) 
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Equation  (77)  is  a  system  of  two  second-order  ordinary  differential  equations. 
Solutions  of  the  form  exp{-jk  z}  are  sought.  Thus,  =  -jk^  an<*  after 
simplification  (77)  becomes 

If  (78)  is  to  have  a  non-trivial  solution,  the  determinant  of  the  coefficient 
matrix  must  be  zero.  This  constraint  yeilds  Fresnel’s  equation  and  permits 
the  solution  for  the  z-comp<5nent  of  the  wave  vector,  k  .  After  some  algebraic 
manipulation,  the  following  quadratic  in  is  obtained. 
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Likewise,  the  tangential  components  of  the  magnetic  field  are  of  the  form: 

-i  4 


H  (z)  = —  7  q.  k.  exp{-jk.z) 

xv  0>U  £  i  v  J  i 


Be  3  4 

H.(z)  =  1  Y.  q.  k  exp{-jk  z} 

111  r 


(83) 


f  ojy 
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The  constants  q.(i  =  1,  2,  3,  4)  can  be  evaluated  by  applying  boundary  condi¬ 
tions  to  the  tangential  components  of  the  electric  and  magnetic  fields. 

Using  (82)  and  (83),  the  tangential  field  components  may  be  expressed 
in  matrix  form  as  follows. 


Q(z)  =  r  e(z)Q 

where 

Q(z)  = 


q2 

q3 

_  q4 
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Hy  (Z) 
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(84) 


! 


Consider  the  layered  shield  depicted  in  Figure  24.  In  each  region  the 
tangential  field 


Region  1 

Region  2 

i 

Region  3 

Figure  24.  Layered  Anisotropic  Shield 

components  are  given  by  (84) 

Region  1.  Q(z)  =  z  <  Z1 

Region  2.  Q(z)  =  T2e2(z)Q2  Zl  <  z  <  Z2  (85) 

Region  3.  Q(z)  =  T3e3(z)Q3  z  >  Z2 

Column  vectors  Q^,  Q  are  chosen  such  that  the  tangential  fields  Q(z)  are 
continuous  throughout  the  three  regions. 

Q(V  =  riGl(Zl)Ql  =  r2e2(Zl)Q2  (86A) 

Q(Z2)  =  r2e2(Z2)Q2  =  r3G3(Z2)Q3  (86B) 

From  (86B) 

Q2  =  [r2e2(Z2))-1Q(Z2)  .  (87) 

Substituting  (87)  into  (86A)  yields 

Q(Zl)  =  r2E2(Zl)[r2e2(Zl)]"lQ(Z2) 

=  r2[e2(Z1)e~1(Z2)]r21Q(Z2)  (88) 

=  tr2e21(z2  -  z1)r21]Q(z2) 

Equation  (88)  specifies  the  transfer  matrix  (note  that  this  transfer 
matrix  is  for  output  variables  in  terms  of  input  variables)  for  the  region 

<  z  <  Z2.  The  quantity  Z2  -  is  simply  the  thickness  of  the  region.  For 
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a  system  of  N  layers,  each  of  thickness  d . ,  the  overall  transfer  matrix  is 
simply  the  product  of  the  transfer  matrices  for  the  individual  layers. 


-1  -1 

t  -  n  r.  c  (d.)r.1 

j-l  LJ  3  J_ 


From  equation  (89)  the  transfer  matrix  of  a  single  anisotropic  layer  of 
thickness  d  can  be  written  in  the  following  form. 


11 

A12 

B11 

B12 

l21 

A22 

B21 

B22 

11 

Ci2 

Du 

°12 

21 

C22 

D21 

°22 

where 


1  2  ^  2  ^  2  A 

^11  =  ^11  *  - 2^  COS  ^  COS  ^  cosh  0^  +  sin  t  cosh  0^ 

cos  U  — 


1  P  2  *  2  ~  2  ^ 

^22  =  ^22  *  - 2 *  s^n  ^  cosh  0^  +  cos  0  cos  r3  cosh  0, 

cos  uL-  ■ 


A21  =  D12  =  - 2^  sin  T"  COS  ^ 

cos  U 


A 

cosh  0^  - 


cosh  0r 


A12  "  °21  =  COS  *  A21 


1  r  4 1  ni 

cos  <p  - 

)s  U*-  cos 


2  *  *  *  2  ^ 
cos  Q  sinh  0^  +  r)^  cos  0  sin  c,  sinh  0^ 


1  2  A  ^1  A  A  /S  A 

^12  =  ^21  **  - 2^  COS  ^  - - X  s^-n^  “  ^2  COS  ^  ^2  S^n  ^  C0S 

cos  U f-  cos  V 


112  *  ~  *  2  A 

B22  *  - 2 a  - X  sin  5  sinh  0^  +  cos  <P  cos  C  sinh  0^ 

cos  U  cos  V 


58 


"11  2~ 
cos  U 


cos  V  2  1  .2  ,  ,  n 

-  cos  r  sinn  0,  +  -  sm  r  sinh  0„ 

Aw  ^  1  A  A  S  2 

n2  cos  <j> 


C12  C21 


2* 
cos  U 


sin  c,  cos  C 


PAC  V  A  1  A  A 

— r —  sinh  0. - cos  <p  sinh  0O 

i  ^  Z 


n. 


n* 


22  2- 
cos  U 


cos  V  .  2  .  “  .  1  3  2  2  r  .  ,  £ 

-  sin  £  sinh  0,  +  —  cos  <p  cos  C  sinh  0O 

12  2 


nn 


n. 


and 


2  "  -  or  .  2 

cos  (p  =  1  -  sin  (p 

e2  o 


2  £  _  or  .2  2 

cos  U  =  1  -  -  sm  d>  cos  C 

e2  ° 


2  ,A.  2  £  or  .  2  ,  .2 

cos  V  =  cos  U  -  -  sin  cp  sm  C 

o 


ei  =  jkid 


02  =  3k2d  ' 


Scattering  parameters  and  hence  reflected  and  transmitted  fields  can 
now  be  calculated  using  Equation  71  of  Section  5.5  of  this  report.  Equations 
(90)  reduce  to  Equation  (66)  of  Section  5.3  when  the  angle  of  incidence,  $  9 
is  set  equal  to  zero. 

6.0  Summary  and  Conclusions 


"ABCD"  and  Scattering1’  matrix  analysis  techniques  for  determining 
shielding  effectiveness  of  isotropic  and  anisotropic  multilayered  shields  are 
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presented.  These  techniques  are  well-suited  to  computer  implementation.  Such 
programs  have  been  generated  and  reported  on  separately. 

Given  appropriate  intrinsic  material  properties  as  input,  the  computer 
simulations  yield  surface  transfer  impedance  and  shielding  effectiveness 
numbers  in  excellent  agreement  with  Boeing  [6]  experimental  data.  The  programs 
can  also  be  used  to  infer  intrinsic  material  parameters  given  measured  surface 
transfer  impedance  and/or  shielding  effectiveness  data. 

From  reported  measurements  on  advanced  composite  materials,  it  is 
apparent  that  the  major  shielding  problem  associated  with  these  materials 
arises  from  seams  and  joints.  Even  in  the  laboratory  under  well  controlled 
conditions,  it  has  proven  exceptionally  difficult  to  obtain  reliable,  repeat- 
able  joints.  Many  of  the  measured  results  appear  erratic  (especially  in 
comparing  between  different  measurement  schemes)  largely  due  to  uncontrollable 
seam  and  joint  leakage  where  the  composite  sample  is  mated  to  the  measurement 
system. 
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APPENDIX  A 


DERIVATION  OF  EQUATION  (9) 

Moser  [7]  and  Bannister  [8,  9]  have  provided  an  integral  equation 
solution  for  the  shielding  effectiveness  of  the  two-loop/ inf inite  flat  plate 
geometry  assuming  uniform  current  in  the  loops.  It  is  assumed  that  the  shield 
can  be  neglected.  The  complete  expression  for  the  shielding  effectiveness  as 
given  by  Bannister  is  shown  in  equation  (A-l) . 
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*  skin  depth  in  the  shield 


U  =  relative  permeability  of  the  shield 

t  =  shield  thickness 
a  =  shield  conductivity 

z  =  r]  r 2  =  center-to~center  separation  of  the  two  Loops 

J^(Aa)  =  Bessel  function  of  order  one  and  argument  (Aa), 
a  =  loop  radius 

A  =  dummy  variable  of  integration 


r  =  J&2  + 


A  .  *  wavelength  in  air 
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Simplified  approximate  forms  of  Moser’s  formula  can  be  derived  under 
certain  conditions  [3,  4],  Let  r1  =  [a^  +  (z  -  t)^]^^  be  termed  the 

measurement  distance.  Then  for  conditions  such  that  r*  <  —  ,  t  t  >  2, 
xrr'  20  r 

Trr*  >  10  and  — —  >  10  are  satisfied,  the  shielding  equation  takes  on  the 
r 

much  simpler  form  given  below. 
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as  is  usually  the  case,  z>>t,  then: 
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If  in  addition  z>>a,  then: 
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MISSION 

of 

Rome  Air  Development  Center 

RAVC  plant  and  execute*  research,  development,  test  and 
detected  acquisition  programs  in  support  of  Command,  Control 
Communications  and  Intelligence.  \C3l)  activities.  Technical 
and  engineering  support  uiithin  areas  of,  technical  competence 
as  provided  to  ESP  Program  Offices  ( POs )  and  other  ESP 
elements.  The  principal  technical  mission  areas  are 
communications,  electromagnetic  guidance  and  control,  sur¬ 
veillance  of  ground  and  aerospace  objects,  intelligence  data 
collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 


